Quark induced excitations of the instanton liquid 



S.V. Molodtsov, A.M. Snigirev''", G.M. Zinovjev^ 

State Research Center, Institute of Theoretical and Experimental Physics, 117259, Moscow, Russia 
^ Skobeltsyn Institute of Nuclear Physics, Moscow State University, 119899, Moscow, Russia 
^Fakultdt fiir Physik, Universitdt Bielefeld, D-33501, Bielefeld, Germany 
^Bogolyubov Institute for Theoretical Physics, 
National Academy of Sciences of Ukraine, UA-03143, Kiev, Ukraine 



O 

Q ' Abstract 

(N : 

^j^, The selfconsistent approach to the quark interactions in the instanton hquid is developed 

^ I within the tadpole approximation calculating the basic functional integral. The effective La- 

' grangian obtained includes colourless scalar field interacting with quarks. The origin of this 

(N 



dynamical field as an interaction carrier in soft momentum region is discussed. 
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g I. Introduction 

'"^ Nowadays, there are no doubts the model of QCD vacuum as the instanton liquid (IL) is the most 
^ ^ practical instrument on the chiral scale of QCD. It provides, as the lattice calculations recently con- 
Q firmed, not only the theoretical background for describing spontaneous chiral symmetry breaking 
(SCSB) but is mostly powerful in the phenomenology of the QCD vacuum and in the physics of light 
quarks while considered to propagate by zero modes arising from instantons. The origin of gluon and 
^ . chiral condensates turns out in this picture easily understandable and both are quantitatively calcu- 
lated getting very realistic values defined by Aqcd and parameters of instanton and anti-instanton 
ensemble, for example, —i (V'V) ~ —(250 MeV)^. Moreover, the scale for dynamical quark masses, 
M ~ 350 MeV, naturally appears and pion decay constant, ~ 100 MeV, is then transparently 
calculated. 

Another significant advantage of this approach is that the initial formulation starts basically from 
the first principles and subsequent approximations being well grounded and reliably controlled are 
plugged in [|I|, 0. It becomes clear especially in recent years when the impressive progress has been 
reached in understanding the instanton physics on the lattice . Further we are summarizing several 
things we have learned thinking of the IL theory and trying to answer the challenging questions. 

Let us start on that stage of the IL approach when its generating functional has already been 
factorized as 

Z = Zg • Z^ , 

where eventually the factor Zg provides nontrivial gluon condensate while the fermion part Z^ is 
responsible to describe the chiral condensate in instanton medium and its excitations. It is usually 
supposed the functional integral of Zg is saturated by the superposition of the pseudo-particle (PP) 
fields which are the Euclidean solutions of the Yang-Mills equations called the (anti-) instantons 



TV 
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Here A^{x; ■ji) denotes the field of a single (anti-)instanton in singular gauge with iNc (for the SU {Nc) 
group) coordinates 7 = (p, z, U) of size p centred at the coordinate z and colour orientation defined 
by the matrix U. The nontrivial bloc of corresponding Nc x A^^, matrices of PP is a part of potential 
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= ^^^^, — 2U'^TaU, y = X - z , a = 1,2,3 , (2) 
9 y y +P 

where are the Pauli matrices, 77 is the 't Hooft symbol 0], (7 is the coupling constant and for 
anti-instanton f] ^ rj. For the sake of simplicity we do not introduce the distinct symbols for 
instanton (A^+) and anti-instanton (A^-) and consider topologically neutral IL with = = N/2. 
Formulating the variational principle the practical estimate of Zg was found p| 

Zg ^ e-<^> 

with the action of IL defined by the following additive functional Q 

{S) = dz dp n{p) s{p) . (3) 



The integration should be performed over the IL volume V along with averaging the action per one 
instanton 

s{p) = Pip) + 5 In(Ap) - In^^JV^ + /^^V' / dpi n{p{)pi, (4) 
weighted with instanton size distribution function 

n{p) = C e"^(^) = C p-5/^27Vcg-/3(p)-.pV^ ^ (5) 

6-4 , II Nc-2 Nf 
b = , 
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where = J dp p n{p)/n = (^ ^^^^ J 5 = I dp n{p) = N/V and Nj is the number of flavours. 
The constant C is defined by the variational maximum principle in the selfconsistent way and 
P{p) = = — IiiCat^ — 61n(Ap) (A = Ap^ = 0.92Ap,y.) with constant Cj^^ depending on the 

renormalization scheme, in particular, here Cj^^ ^ n'^^N ^"^^^^ 1 ^j^^^2) ! ' '^^^ parameters (3 = l3{p) 

and (3 = (3 + In Cn^ are fixed at the characteristic scale p (an average instanton size). The constant 
27 A^ 

= j^'i ^ ^ T^"^ characterizes, in a sense, the PP interaction and Eqs. (^,(§) and (|^) describe the 

equilibrium state of IL. The minor modification of variational maximum principle (see Appendix) 

leads to the explicit form of the mean instanton size pA = exp j^^^^jj a-^^d, therefore, to the 

direct definition of the IL parameters unlike the conventional variational principle |^ which allows 
one to extract those parameters solving numerically the transcendental equation only. 

The quark fields are considered to be influenced by the certain stochastic ensemble of PPs, Eq. 
(|I|), while calculating the quark determinant 



Z.^ c:i ( Dtfj^D^ {{ e^^^'^'^^^ )) 
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^ In fact, the additive property results from the supposed homogeneity of vacuum wave function in metric space. 
Eq. looks like a formula of classical physics although it describes the ground state of quantum instanton ensemble. 
Intuitively clear, this definition will be still valid even when the wave function is nonhomogeneous with the nonuni- 
formity scale essentially exceeding average instanton size or, precisely speaking, being larger (or of the order) than 
average size of characteristic saturating field configuration. Then each instanton liquid element of such a distinctive 
size will provide a partial contribution depending on the current state of IL, see next Section. 



Besides, dealing with the dilute IL (small characteristic packing fraction parameter n p^) one neglects 
the correlations between PPs and utilizes the approximation of iVg — >■ oo where the planar diagrams 
only survive. In addition the fermion field action is approached by the zero modes which are the 
solutions of the Dirac equation i{]D{A±) + m) $± = in the field of (anti-)instanton A±, i.e. 



P 



V^7r|x|(x2+p2)3/2 
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with the colour c,d and the Lorentz i,j indices and the antisymmetric tensor e. In particular, at 
Nf = 1 the quark determinant reads 



= i dz dU dp n{p)/n / dxdy ip\x) idx^±{x — z) — z) idy ip{y) , 



where the factor R makes the result dimensionless and is also fixed by the saddle point calculation. 
Much more accurate result for the Green function of quark in the PP ensemble allows one to 
calculate the generating functional even beyond the chiral limit and the simple zero mode approxi- 
mation turns out amazingly fruitful again to develop (quantitatively) the low energy phenomenology 
of light quarks . 

Thus, the IL approach at the present stage of its development looks very indicative, well the- 
oretically grounded and reasonably adjusted phenomenologically. The proper form of generating 
functional obtained and its reasonable parameter dependence provide with enough predictive power 
and justify, hence, the approximations made. Generally, it leads to conclude the quark feedback 
upon the instanton background is pretty limp and could be perturbatively incorporated as a small 
variation of instanton liquid parameters 5n and 5p around their equilibrium values of n and p being 
in full analogy with the description of chiral condensate excitations. Indeed, the result of nontrivial 
calculation of the functional integral (treating the zero quark modes in the fermion determinant 
substantially from physical view point) takes down to 'encoding' the IL state with just those two 
parameters. Moreover, the IL density appears in the approach via the packing fraction parameter 
np^ only (clear from dimensional analysis) what means one independent parameter existing in prac- 
tice. It is just the instanton size. The analysis of the quark and IL interaction is addressed in this 
paper developing our idea ||^ of phononlike excitations of IL resulting from the adiabatic changes 
of the instanton size. Hence, we describe further the quark feedback dealing with these deformable 
(anti-)instanton configurations which are the field configurations Eq. (Q) characterized by the size p 
depending on x and z, i.e. p p{x, z). 

The paper is organized as follows: in Section II we discuss the modification of quark determinant 
when the functional integral is saturated by the deformable modes at one flavour. Then in Section 
III we develop the proper calculation (tadpole approximation) which is based on the saddle point 
method. Section IV is devoted to the generalization for the multifiavour picture. The paper includes 
also Appendix where the fault finding reader gets a chance to control the explicit formulae for the 
IL parameters and to improve our calculations if is able. 



II. Dynamics of phononlike excitations 

Apparently, the essence of what we discuss here could be illuminated in the following way. Saddle 
point method calculation of the functional integral implies the treatment of the action extremals 
which are the solutions of classical field equations. For the case in hands the action S[A,%p\ip\ is 
constructed including the gluon fields A^, (anti-)quarks ip\ip and extremals which are given by the 
solutions of the consistent system of the Yang-Mills and Dirac equations. As a trial configuration 



in the IL theory the superposition of (anti-)instantons which is the approximate solution of the 
Yang-Mills equations (with no backreaction of the quark fields) and a background field for the 
Dirac equation is simultaneously considered. We believe it is reasonable to make use the deformable 
(crumpled) (anti-)instantons A±{x; 7(0:)) as the saturating configurations. They just admit of varying 
the parameters 7(x) of the Yang-Mills sector of the initial consistent system in order to describe the 
influence of quark fields in the appropriate variables for the quark determinant. 

Taking the action in the form S[A±{x,'y{x)),'^|J\^p] we would receive the corresponding varia- 
tional equation for the deformation field 7(2;) which approaches most optimally (as to the action 
extremum) PP at nonzero quark fields. In the field theory, for example, the monopole scattering |^ , 
the Abrikosov vortices scattering |]I0| are treated in a similar way. Actually, the unstationary picture 



seems more adequate as a specific role of instanton-anti-instanton annihilation channel prompts and 
an investigation of two-body colour problem teaches |TT . 



In practice, for IL we avoid the difficulties which come with solving the variational equations 
if we consider the long wavelength excitations only (with the wave length A much larger than the 
characteristic instanton size p). Indeed, it can be done because we are searching the kinetic energy 
of the deformation fields | (the one particle contributions) and consider the pair interaction which 
develops the contact interaction form being calculated in the adiabatic regime [Q . 

Let us remember here that deriving Eq. (|^) we should average over the instanton positions in 
a metric space. Clearly, the characteristic size of the domain L which has to be taken into account 
should exceed the mean instanton size p. But at the same time it should not be too large because 
the far ranged elements of IL are not 'causally' related. The ensemble wave function is expected 
to be homogeneous (every PP contributes to the functional integral being weighted with a factor 
proportional to ~ 1/V, V = L^) on this scale. The characteristic configuration which saturates the 
functional integral is taken as the superposition Eq. (|l|) with number of PP in the volume V. It 
is easy to understand that because of an additivity of the functional, Eq.(|^) describes properly even 
non-equilibrium states of IL when the distribution function n{p) does not coincide with the vacuum 
one and, moreover, allows us to generalize this action for the non-homogeneous liquid when the size of 
the homogeneity obeys the obvious requirement X > L > p. In particular, the ensemble of deformable 
PPs as saturating the functional integral obeys these constraints if the corrections coming with the 
deformation fields g^iy are small and the long wavelength excitations (on the scale of the instanton size 
dp{x, z) 



p, I.e. 



dx 



essential simplification to define the field in the center of instanton ^V. ' — ^ ~ ' — '- 
^ ox ox 



^ 0(1)) of initial instanton fields Gf^,y are mainly taken into consideration. Such 
a condition of smoothness or the adiabatic change of instanton size dictates, in practice, another 

as a 

ox ox 
characteristic deformation field in further calculations. 

In principle, the correction to the instanton field produced by the quark presence might be 
calculated in its most general form, if the gluon Green function in the instanton medium is known, 
as 

a';{x, z) = jdi D;1{x -z,^-z) Jti^ - z^) , (8) 

where is the current of external quark source, belongs to the region of long wavelength distur- 
bance and, at last, is the Green function of PP in the instanton medium. In fact, this Green 
function (even when considered in the field of one instanton) is not well defined |^ but, seems, for 
the case in hands we could develop the selfconsistent way to calculate the regularized Green function. 
The nonsingular propagator behaviour in the soft momentum region is defined by the mass gap of 
phononlike excitations. Fortunately, the exact form of the Green function occurs to be not of great 
importance here (we are planning to return to the problem of regularized Green function calculation 
in the forthcoming publication). In the coordinate space it is peaked around the average PP size 

^Then we are allowed to take the slowly changing deformation field beyond the integral while calculating the action 
of deformed instanton. 



being in nonperturbative regime and, hence, integral Eq. (|^) is estimated to be 



a;ix,z)c^D';tix-z) J',iz-z^) . (9) 
On the other hand, the exact instanton definition in the singular gauge, 



g dxy 

leads to the following correction to the instanton potential 




<(^) = ^;.(^,^)^, (10) 



where H°;„(x,z) = — a^^^^'/^^ — ^. Confining within the precision accepted here we could take 

p{x, z) ~ p, i.e. H^X^, z) = H^^{x - z). 

If we compare now both definitions of the correction we are capable to get immediately for 
the following equation 



- z) = Df^{x - z) Jl{z - z^) . 



In fact, the current J might be taken constant in the long wavelength approximation. Then neglecting 
the current gradients one is allowed to change the derivatives to obtain the following estimate of the 

deformation velocity 



dp{x, z) 




dp{z) 


dx 




dz 



which looks to be well justified since there is no other 

fields in the problem at all (in the adiabatic approximation). The contribution of deformed (anti- 
)instantons to the functional integral (when the corrections coming from the PP deformation fields 
are absorbed) may be estimated as 

(5) f dz f dp n{p) 




where k is the kinetic coefficient being derived within the quasiclassical approach. Our estimate of 
it gives the value of a few single instanton actions k ^ c j3 with the coefficient c ~ 1.5 ^ 6 depending 
quantitatively on the ansatz supposed for the saturating configurations. Although this estimate is 
not much meaningful because there is no the vital k dependence eventually (becomes shortly clear) 
and the kinetic term could be introduced phenomenologically. Thus, this coefficient should be fixed 
on a characteristic scale, for example k ~ ^(p), if we not are planning to be beyond the precision 
peculiar to the approach. Actually, it means adding the small contribution of kinetic energy type to 
the action per one instanton only. Such a term results from the scalar field of deformations and affects 
negligibly the pre-exponential factors of the functional integral. In one's turn the pre-exponential 
factors do the negligible influence on the kinetic term as well. The deformation flelds induced by 
the dilations and rotations in the isotopic space result in the singular and trivial kinetic coefficients, 
respectively §]. Then to evaluate the mass scale related to those modes one needs the heuristic ideas 
rooted, apparently, beyond the conventional IL and SCSB. 

If we strive for being within the approximation developed we should retain the small terms of 

dsjp) 



the second order in deviation from the point of action minimum 



dp 



only expecting the 



approximate validity of 



P=Pc 



s{p) ^ s{p) + , (11) 



where s^'^\p) 



<fs{p) 



Pc 



= and the scalar field 
p2 



deviations from the equilibrium value of pc = P ( ■'- ~ 2z7 
field is described by the following Lagrangian density 



Sp = p — Pc — p — p is the field of 
p. Consequently, the deformation 



C 




2,„2 



with the mass gap of the phononlike excitations 

= — = 

which is estimated for IL with = 3, for example, in the quenched approximation to be 

M ^ 1.21 A 

if c = 4, p A ^ 0.37, (3 ^ 17.5, n A~^ ^ 0.44 (for the details see the tables of Appendix). 
Changing the variables we obtain the gluon part of the generating functional as 



4~ 



SA 



6ip ■ 



UK 

exp < — — I dz 



'dip' 
dz 



with the Jacobian 



5A 

dip ' 



corresponding those new variables introduced. Let us notice that the 
Decause of incomplete set of the transformations shown. However, in the 



latter looks pretty formal 
adiabatic approximation, as was mentioned above, the preexponentional Jacobian contribution to 
the action f\ being a c-number should be omitted. 

Analyzing the modifications which arise now in the quark determinant we take into account 
the variation of fermion zero modes resulting from the instanton size perturbed 



$±(X - 2,p) ~ - Z,Pc) + ^±\x - Z,Pc) Sp{x,z) 



where $£^('U, p^ 



d<^±{u,p) 
9p 



and because of the adiabaticity it is valid 6p{x, z) ~ (5p(z, z) 



P=Pc 



ip{z). The additional contributions of scalar field generate the corresponding corrections in the 
factors of the kernels of Eq. (0) which are treated in the linear approximation in ip and taking 
approximately pc = p, i.e. 



id^^±{x - z,p) ^\{y - z,p) idy ^ T±{x,y,z,p) + r2^(x, p) ip{z) , 



'121 



here we introduced the notations 



r±(a;,y,z,p) = id^^±{x - z,p)^l^{y - z,p){-idy) 



yj ' 



r'l\x,y,z,p)=tdM'ix-z,p)<^>l{y-z,p) {-tdy)+td,<l>±{x-z,p) <l>l"{y-z,p) {-tdy) 



.(1) 



t{i)/ 



■^Generally, the deformation field p^, and integration variable a° (^0|) are related via the rotation matrix: a° = 
^ab^'liiyix — z)p^ and in the long-length wave approximation $ might be constant $J|^i^(0) {x ~ z). With the rotation 
matrix spanning the colour field = ri"^a^ on the fixed axis we can conclude that the vectors and are, in 
fact, in one to one correspondence (of course, being within one loop approximation and up to this unessential colour 
rotation). Thus, the Jacobian contribution turns out to be an unessential constant. 



with —idy left acting operator (the gradients of scalar field (p are negligible according to the adia- 
baticity assumption again). It is a simple matter to verify that the right hand side of Eq. ([I2|) , being 
integrated over dzdU, generates the following kernel (in the momentum space) 



(13) 



with 7o(fc,A;) = G^{k) , G{k) = 2TTpF{kp/2) , -fi{k,l) = G{k)G'{l) + G'{k)G{l) 



G'{k) 



dG{k) 
dp 



F{x) = 2x [Io{x)Ki{x) - Ii{x)Koix)] - 2 Ii{x)Ki{x), where h, Ki {i = 0, 1) 



p=p 



are the modified Bessel functions. 

The functional integral of Eq. (|^) including the phononlike component may be exponentiated in 
the momentum space with the auxiliary integration over the A-parameter (see, for example 0) 



X / Dijj'^Dip exp 



iXVR 
dkdl 



{2nf 



iX 



iX 



{27:)%k -l)\-k + — 7o(fc, k) + — 7i(fc, /) ^{k - I) 



N. 



(we dropped out the factor normalizing to the free Lagrangian everywhere). It is pertinent to mention 
here the Diakonov-Petrov result comes to the play precisely if the scalar field is switched off. 

In order to avoid a lot of the needless coefficients in the further formulae we introduce the 
dimensionless variables (momenta, masses and vertices) 



kp 
~2 



Mp 



M , 7o p^7o 



{np K 



(14) 



the fields in turn 



ip{k) {nK)-^l^p^^{k) , i){k) p^'^i){k) , (15) 
Xtf' 

and eventually for A we are using p = -cjjpf - Then the generating functional takes the following form 

dk 1 



Z ~ J dp Zg J Dip"^ Dip Dip exp |np^ |^ln 
dkdl 
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np 
p 



vr^ 2 



ip{-k) 4 [k' + M'] ip{k) \ X 



(16) 



X exp 



TT 



ij\k)2 n^6{k-l){-k + ipjo{k,k))+ip-fi{k,l) ^{k-l) ij{l) 



where Zg is a part of gluon component of the generating functional which survives after expanding 
the action per one instanton Eq. (pH]). The functional obtained describes the IL state infiuenced 
by the quarks when all the terms containing the scalar field are collected (see also Appendix). As 
mentioned above we believe this infiuence analogous to the back impact of phononlike deformations 
on the quark determinant does not considerably change the numerical results of the IL and SCSB 



^In the metric space we have the nonlocal Lagrangian of the phononlike deformations (p{z) interacting with the 
quark fields ip'^ , ip, i.e. 



dx ip'^ {x)idxip{x) 



dz 



dip 
dz 



+ M^ip'{z) } + 



dxdydz dU 'iJj^x){T±{x,y, z, p) + T''l\x,y, z, p) ip{z)}'ip{y) 



The physical meaning of the basic phenomenon behind this Lagrangian seems pretty transparent. The propagation 
of quark fields through the instanton medium is accompanied by the IL disturbance (the analogy with well known 
polaron problem embarrasses us strongly in this point). 



theory. Noninteracting part of phononlike excitation Lagrangian characterizes the IL reaction on 
the external long-length wave perturbation and, apparently, is its general feature independent of 
perturbating field nature. 

The presence of quark condensate makes hint for the appropriate scheme of approximate calcu- 
lation of the generating functional 

III. Tadpole approximation 

The formal integration over the scalar field leads us to the four fermion interaction and the functional 
integral can not be calculated exactly. However, due to smallness of scalar field corrections we may 
find the effective Lagrangian substituting the condensate value in lieu of one of the pairs of quark 
lines (see Fig. la.) 

i'^ik)ij{l) {i)\k)i){l)) = -T^H{k - I) Tr S{k) 

where S{k) is the quark Green function. In such an approach the diagram with four fermion lines in 
the lowest order of the perturbation theory in /x is reduced to the two-legs diagram with one tadpole 
contribution (there are two such contributions because of two possible ways of pairing) 

r rlhrll rlh'ril' 

2 ^,{k,l) 7i(A;',n ^I^Hk)m ^Hk'm') - I) <p{k' - I')) ^ 

^ 4 /.^ / 4 k) ^P\k)m I 4 7i(^ Tt S{1) DiO) , 
where the natural pairing definition was introduced 

{<f{k) ^{l)) = 7r'5{k + I) Dik), D{k) = ^^i^2\m2^ ■ 

It is obvious the factors surrounding if)^ {k)'il){k) have a meaning of an additional contribution to the 
dynamical mass 

m{k) = 7i(^, k) {-2if,) / - I) Tr S{1) D{0) , (17) 

J 77 

(the initial mass term contains the factor 2 when the dimensionless variables are utilized, i.e. takes 
a form 2im.) 

The contribution of the graph with all the quark lines paired (see Fig. lb) 



-2 



/ ^ 7i(^, k) Tr S{k) 7r%0)D{0) = ^ - f^^,{k,k)Tr S{k) 
J 7T I p V J 7: 



p^ V K 



dk 



2 1/ 

together with contribution of the graph (Fig. Ic) 

r dkdl 

2 /iM ^ Tr7i(A;, /) 7i(/, k) S{k) S{1) D{k - I) , 

J TT 

should be taken into account at the same order of the p expansion while calculating the saddle point 
equation. Here we used the natural regularization of the 5-function S{0) — in the dimensionless 
units. Then the quark determinant after integrating over the scalar field reads 



Z ~ [dp [ Di/j^Di/j exp \np^ ( In ^ - iW ^ ^ p'' c^ip)- 
J J \ \ p I np u p 



1)+— 4^-/" C^(/^)- 



- ^^'''V / — ^?<*- " (t^ I rV-))(f '+ r-(0) ^''^ - " + ^? /^^ ^ iuM + .rwi| . 

where the vertex function is defined as 

T{k) = /i 7o(fc, fc) +m(fc) , 
and we introduced the function c(/i) convenient for the practical calculations 

i{np^Ky/'^ f dk 



As seen from Eq. (|18D the Green function of the quark field is selfconsistently defined by the following 
equation 

2 [-k + iV{k)] S{k) = -1 . 

Searching the solution in the form 

S{k) = A{k) k + i B{k) , 

we get 

1 1 1 



Mk) = ^ T^-FoTTT , B{k) = - 



2 k^ + T\k) ' ' ' 2 k^ + T\k) ■ 

Using Eq. (|l^) and the definitions of T{k) and B{k) we have the complete integral equation 

r{k) = p 'joik, k) + AT,- fc) / 4 7i(/, W^T^Tl- 



which drives to have the convenient representation of the solution 

r(A;) = II 7o(A;, k) + ^^^777^- /i^ c(/i) 7i(A;, fc) . 
What concerns the function c(/i) it is not a great deal to obtain 

and, therefore, the complete integral equation for the function c(/i) which is shown in Fig. 2 for 
N j = 1. Let us underline the A^j-dependence of the c(/i) function in the interval of fi determined 
by saddle point value is unessential. Then we easily obtain for the additional contribution to the 
dynamical quark mass 

and see the cancellation of the kinetic coefficient k in m if we remember the definition (P^). Thus, 
it means the precise value of the coefficient is unessential as declared. 

We have the following equation for the saddle point of the functional of Eq. 



p 'J e + T\k) ' V 71^ \ {e + T\k)){P + T\l)) 



(20) 



where the prime is attributed to the differentiation in fi. Being interested in receiving a closed 
equation we need to know the derivative c'(/x) too. Then the definition of c(/x) above allows us to 
have 

(1 - fx^A{fx)) c'(/i) = 2/i A{fx) c{fi) + B{fx) , 



where 



B 



N,K f dk -ffjk) e - r^jk) 
e + r^{k) e + r\k) ' 

2{np'Hy/^ f dk 7i(A-) r'{k) 



Ai„) = ail.) ^ J 



2 ' 



and 



a{ii) — 1 — /J, 



2iV, r(z/+l/2) c(/x) 



f3e v^l^Viy) 



np I np 



Thereby the saddle point equation absorbs the effect of the IL parameter modification {n{ii)) pro- 
duced by equilibrium instanton size shift Pc P which comes in the leading order from simple tadpole 
graph 



^ 7i(A;, (-TT^) 5{k - I) Tr S{k) ip{k - I) = A ■ ip{0) , 

TT 

-fi{k,k) Tr S{k) 

TT 



(21) 



{np k) 



(remind here (p — p — Pc, and </?(0) — J dz (f{z) is the scalar field in momentum representation). 

In the Tabic 1 the numerical results (M.S.Z.) are shown for Nf — 1 comparing to those Diakonov 
and Petrov (D.P.) 

Table 1. 



D.P. 


M.S.Z. 


// M(0) -iiijj^) 


p M(0) -i(V'V) 


5.27-10-=^ 359 -(333)=^ 


4.81-10-3 362 -(322)3 



The parameters indicated in the Table 1 are 
the dynamical quark mass 

M(0) = 2r(o) 

and the quark condensate 



1 



[MeV] 



dk T{k) 



1 



[MeV] 3 . 



Through this paper the value of rcnormalization constant is fixed by A = 280 MeV. Then the IL 
parameters are slightly different from their conventional values p ~ (600 MeV)~^, R ~ (200 MeV)~^ 
(see the corresponding tables in Appendix). However, with the minor A variation the parameters 
could be optimally fitted. As expected the change of quark condensate is insignificant, the order of 
several MeV, what hints the existence of new soft energy scale estabhshed by the disturbance which 
accompanies the quark propagation through the instanton medium. 



IV. Multiflavour approach 



In order to match the approach developed to phenomenological estimates we need the generahzation 
for Nf > 1. Then the quark determinant becomes 



J f=i 



' Y+ 
VR^f 



Nf 



N+ 



Y- 



N- 



VR^f 



= i f dz dU dp n{p)/n Y[ / dxfdyf ^f{xf) idx^.^±{xf - z) $±(?// - z) idy^ i^fiVf) ■ 

/=i 

With phononhke component included every pair of the zero modes ~ $ <|)t acquires the additional 
term similar to Eq. ( [T^ ) . The appropriate transformation driving the factors F ^ to their determinant 
forms is still valid here since the correction term differs from the basic one with the scalar field 
if. The complete integration over dz leads (in the adiabatic approximation ip{x,z) — * <^{z)) to the 
transparent Lagrangian form with the momentum conservation of all interacting particles. Besides, 
we keep the main terms of Y^ expanding in (p. The quark zero modes generate the factor similar to 

Eq. ( p!3D with being changed by the factor ( j and then in the leading Nc order we have 



Y-=[-) jdz^t{^J\z) 



(2vr) 



dp 



While providing the functional with the Gaussian form we perform the integration over the auxiliary 
parameter A together with the bosonization resulting in the integration over the auxiliary matrix 
Nf X Nf meson fields [|12| 



exp 



, , , ^ J 
Adetl- 



dM exp { I Tr[MJ] - {Nf - 1) 



'Aet[MN^ 



As a result the generating functional may be written as 
Z = J^Z'; exp(-iVlnA) J exp|-| 

■ j DMl,r exp J dz |-(iV/ - 1) 



dk UK 



(27r)^ 2 

deilMLNcW^ , (det[MRN^]\~- 



+ 



A 



(22) 



■ j D^^Dtfj expl^J ^^Y.^\[k){-k)^^f{k) + i j dz{Tr[MLJ^]+Tr[MRj- 

Now scalar field interacts with the quarks of the different flavours, nevertheless, the dominant con- 
tribution is expected from the tadpole graphs where any pair of the quark fields is taken in the 
condensate approximation as done at Nf = 1 

^}{k)^g{l) - (^}(fc)^,(/)) = -7r%g 6{k - I) Tr S{k) . 

As for the condensate itself we obtain it as the nontrivial solution of saddle point equation. For 
example, it is 

iML,R)fg = M 6fg 



for the diagonal meson fields. The dimensionless convenient variables (in addition to Eqs. (0), (|T5|)) 
are the following 



M 



P 



(2iVeP) 



Nf 



9 ■ 



Then the effective action {Z = J dgdfi exp{— Vg//}) in new designations has the form 



y^ff = N{Nf - 1) In^? - -^{Nf - 1)^^ 



+2N,Nffi- 4 
P 



p np V 



V2NjNl , 

— -4 P C in) + 



TT" - ■(fc2 + r2(A;))(/2 + r2(/)) 

and the saddle point equation reads 



D{k -I)- 2NjN, _4 



V r dk 



71 



(23) 

\n{e + T\k)} , 



np 
P 



dkdl \ ii^-il{k,l)[{kl) -V{k)V{ 



D{k - I) - {np%\n 



TT 



(P + r2(A;))(/2 + r2(/)) 



np 



7.4 \ 1/^/ 



(24) 



. 



The additional contribution to the dynamical quark mass in Eq. (p!9| ) gains the factor Nf because 
the scalar nature of the phononlike field requires to match the tadpole quark field condensates of all 
Nf flavours to each vertex 

[np^Kyi^ V 

Table 2 complements the Table 1 with the calculations at N f = 2 

Table 2. 



D.P. 


M.S.Z. 


/i M(0) -i(V'V) K 


M(o) f; 


4.83-10-3 386 -(381)3 ^35 m 


3.26-10-3 298 -(335)3 ^08 90 



where [MeV] is the pion decay constant and 

N^Nf rdkr\k)~^j^\k)m + ^i^'ik)r (I 



TC 



[MeV] is its approximated form 



y2 



N,Nf 



{e + T\k)f 



dk T\k) 

l^ik^ + r^k))' 



1 

.p. 



dT{k) 



here T'(k) = ^^Jj^^ ? and the condensate —i{^^ip) is implied for the quarks of each flavour. 

The light particle introduced and which imitates the scalar glueball properties does not affect 
significantly the SCSB parameters and correctly describes the soft pion excitations of quark conden- 
sate. Meanwhile, the experimental status of this light scalar glueball is very vague. We believe the 
phononlike excitations could manifest themselves being mixed, for instance, with the excitations of 
the quark condensate in the scalar channel. 



VI. Conclusion 



In this paper we have developed the consistent approach to describe the interaction of quarks with 
IL. Theoretically, it is based (and justified) on the particular choice of the configurations saturating 
the functional integral what is not merely a technical exercise. They are the deformable (crumpled) 
(anti-)instantons with the variable parameters 7(x) and in the concrete treatment of this paper we 
play with the variation of the PP size p{x,z). In a sense, such an ansatz is strongly motivated by 
the form of quark determinant which is solely dependent on the average instanton size in the SCSB 
theory. We have demonstrated that in the long-length wave approximation the variational problem 
of the deformation field optimization turns into the construction of effective Lagrangian for the scalar 
phononlike if and quark fields with the Yukawa interaction. Physically, it allows us to analyze the 
backreaction of quarks on the instanton vacuum. We have pointed out this infiuence on the IL 
parameters as negligible. The modification of the SCSB parameters turns out pretty poor as well. In 
particular, the scale of quark condensate change amounts to a few MeV only. Nevertheless, switching 
on the phononlike excitations of IL leads to several qualitatively new and interesting effects. The 
propagation of the quark condensate disturbances over IL happens in this approach to be in close 
analogy with well-known polaron problem. We imply a necessity to take into account the medium 
feedback while elementary excitations propagating. Besides, it hints that fitting the parameters p, n 
and renormalization constant A all together with the alteration of s{p) profile function we might 
achieve suitable agreement not only in the order of magnitude. The difficulties which confronted us 
here illuminate the fundamental problem of gluon field penetration into the vacuum (the instanton 
vacuum in this particular case) as the most principle one. Indeed, it is a real challenge to answer 
the question about the strong interaction carrier in the soft momentum region. Perhaps, the light 
particle of scalar glueball properties which appears inherently in our approach and should manifest 
itself in the mixture with the excitation of quark condensate in scalar channel (cr-meson) is not bad 
candidate for that role. By the way, it could be experimentally observed as a wide resonance. 

Summarizing, we understand our calculation can not pretend to the precise quantitative agree- 
ment with experimental data and see many things to be done. We are planning shortly to consider 
the problem of instanton profile [O, to make more realistic description of the PP interaction and 



to push our ansatz beyond the long-length wave approximation analyzing more precisely 'instanton 
6 A 



Jacobian' -r— 

dip 
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Appendix 

The contribution of the quark determinant to the IL action is given by the tadpole diagram Eq. (pl|) 
which takes the following form when returned to the dimensional variables (see, Eq. (|T5|)) 

nnf/^ _4 .1/2 /■ , n{p) f dz p{z)-pc 



= A (/.(O) = A {np^Ky/^ I dp"^ I , 

p J n J p p 



^ We receive the effective Lagrangian of Yukawa type wliich admits tfie estimate of quark-anti-quark bound state 

while adapted to the non-rclativistic approximation and the inequality 1^} ^ Tl^^*^ > 2 if valid signals its appear- 

47r M np K 

ance. In the problem under consideration the left hand side of the inequality is 0(1). 



Then the IL action, Eq. (j^), acquires the additional term 



{S) = /rfzn|(s)-(A'^)| , 

where A' = /i^c(/i) and the mean action per one instanton is given by the following functional 

(si) = J dp si(p)n{p)/n with 

si(p) = Pip) + 5 In(Ap) - ln^2iv= + V - A' (p - pj/p • 
In order to evaluate the equilibrium parameters of IL we treat the maximum principle 

adapting it to the simplest version (when the approximating functional is trivial 5*0 = 0). In a 
sense, this choice of the approximating functional should be a little 'worse' than in Ref. [0. Its 
only advantage comes from the possibility to get the explicit formulae for the IL parameters in lieu 
of solving the complicated transcendental equation. In equilibrium the instanton size distribution 
function n(p) should be dependent on the IL action only, i.e. n(p) = Ce~^^^^ where C is a certain 
constant. This argument corresponds to the maximum principle of Ref. 0. Indeed, if one is going to 
approach the functional (H) as a local form (s) = J dp si{p)n{p)/n where si(p) = /3(p) + 51n(Ap) — 
InP"^^" + P^'^p'^np'^, it makes the approach self consistent. The functional difference (s) — {si) = 
J dp {s(p) — Si{p)}e~^^pyn being varied over s(p) leads then to the result s(p) = Si(p) + const 
keeping into the mind an arbitrary normalization. The maximum principle results in getting the 
mean action per one instanton as the IL parameters function, for instance, average instanton size p. 
The corrections generated by the 'shifting' terms turn out to be small and we consider them in the 
linear approximation in the deviation A. The following schematic expansion exhibits how the major 
contribution appears 



here 5(A) stands for a certain small 'shifting' contribution and s is the action generated by the gluon 
component only. The last term in Eq. (^) is small comparing to the first one and we ignore it. 
Then it is clear that evaluating the mean action per one instanton is permissible to hold the gluon 
condensate contribution s only (without the 'shifting' term 6) in the exponential. Hence, we have for 
the mean action per one instanton (si) = J dp si(p)no(p)/no, and no(p) is the distribution function 
which does not include the 'shifting' term |. It is possible to obtain for the average squared instanton 
size and the IL density that f\ 

2^ f A'r(z/ + 1/2)1 f ,,r(z/+l/2)l , , 



The 'shifting' term changes the mass of phononhke excitation insignificantly. The equilibrium instanton size as 

= is equal to = (a+A' f3) p, a ^ (l - ^y^\ /3 ^ ^ 1^1 - a 



dictated by the condition 

and the second derivative of action in the equilibrium point equals to s (pc) — =^ + =^ ^' { ^^3/2p|^^"^ ^ ^2va 
Another source of corrections to the kinetic coefficient appears while one considers the instanton profile change A — > 

. This mode could appear within the superposition ansatz Eq. 

(|^) and leads to the modification of quark zero mode {D{A + a)ijj — 0). Fortunately, both corrections to the kinetic 
term are numerically small. 

^In order not to overload the formulae with the factors making the results dimensionless, which are proportional 
to the powers of A, we drop them out hoping it does not lead to the misunderstandings. 



(J oi X z 

A + a where the field of correction is a ~ '(Jx — 



n = C CnJ 



2r 



2u ' 



where the parameter equals to 



(27) 
(28) 



P 



Expanding Inp = Inp + ^ ^ ^ + -h^^^—^r^ — h ■ ■ ■ and using Eq. (p6D we can show that 



p ' 5 
/rfp rao(p) Inp 



Inp + $i(z/) 



/rfp no(p)p 



p + <l>2(z/) 



Jdpnoip) ' Jdpnoip) 

where $i, $2 are the certain function of u independent of p. Besides, the average squared instanton 
size within the precision accepted obeys the equahty r^p^ = $(z/), and $(z^) is the function of u only. 
Then the mean action per one instanton looks like 

{si) = -2Nc\n(3 + {2u - 1) \np + F{u) 

(F(i/) is again the function of u only and its explicit form is unessential here). Calculating its 
maximum in p we receive 

2Nr 1 ^ 2bNr 



p = exp 



-InC 



2z/ - IJ ' ' 2z/- 1 
Eqs. ([26|), (|28|) allows us to receive the equation for packing fraction parameter 

u 4 A r(z/ + l/2) 



(np 



np 



(29) 



and for positive root we find 



SNc 



n = u 



rv + V2)l 



and handling Eq. ( pTD we determine the constant C. 

The IL parameters come about close to the parameter values of the Diakonov-Petrov approach 
[0 and are shown in the following Table 

Table 3. 



D.P. 


M.S.Z. 


Nf pA n/A^ (3 


Nf pA n/A^ (3 


0.37 0.44 17.48 

1 0.30 0.81 18.86 

2 0.24 1.59 20.12 


0.37 0.48 17.48 

1 0.33 0.70 18.11 

2 0.28 1.13 18.91 



Here Nf is the number of flavours {Nf = corresponds to the quenched approximation) and N^ = 3. 
It is curious to notice the quark influence on the IL equilibrium state provokes the increase of the IL 
density. 

In Table 4 we demonstrate the mass gap magnitude M and the wave length in the 'temporal' 
direction A4 = M~^. To make it more indicative we show also the average distance between PPs 
from which it is clear, indeed, that the adiabatic approximation A>L~^>pis valid for the 
long-length wave excitations of the vr-meson type. All the parameters are taken at k = 4(3 but the 
primed ones correspond to the kinetic term value of k = 6/3. 

Table 4. 



Nf 




AA 


M'A-i 


A'A 


RA 





1.21 


0.83 


0.99 


1.01 


1.2 


1 


1.34 


0.75 


1.09 


0.91 


1.1 


2 


1.45 


0.69 


1.18 


0.84 


0.97 



here R 



n 



'^1^ is the distance between PPs. 
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Figure 1: The tadpole graphs. 




0^000 0.005 0.010 0.015 0.020 



Figure 2: The function c(/x) ai Nf = 1. 



